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Lorentz invariance violation in top-down scenarios of ultrahigh energy cosmic ray 

creation 
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The violation of Lorentz invariance (LI) has been invoked in a number of ways to explain issues 
dealing with ultrahigh energy cosmic ray (UHECR) production and propagation. These treatments, 
however, have mostly been limited to examples in the proton-neutron system and photon-electron 
system. In this paper we show how a broader violation of Lorentz invariance would allow for a series 
of previously forbidden decays to occur, and how that could lead to UHECR primaries being heavy 
baryonic states or Higgs bosons. 

PACS numbers: 11.30.Cp, 96.40.De, 98.70.Sa 



I. INTRODUCTION 

UHECRs are an enduring mystery. Without the in- 
troduction of new particles or interactions, evading the 
Greisen-Zatsepin-Kuz 'min (GZK) cutoff E, requires 
unidentified nearby sources. Even without the GZK cut- 
off, the "bottom-up" approach faces the challenge of find- 
ing in nature an accelerator capable of energies in excess 
of 10 20 eV. 

"Top-down" scenarios assume that the UHECRs result 
from the fragmentation of a ultra-high energy hadronic 
jet produced by cosmic strings 0] or by the decay of 
a supermassive particle 4J. In the supermassive par- 
ticle (wimpzilla) scenario, the UHECRs are of galactic 
origin, resulting from the decay of relic supermassive 
(M > 10 13 GeV) particles. Wimpzillas can be produced 
copiously in the early universe H , thus solving the energy 
problem. Since they would cluster in the dark matter 
halo of our galaxy 0, 0, y| , they also solve the distance 
problem. Detailed analysis of these decays, however, 
show that at high energy top-down scenarios produce 
more photons than protons || fToL [HI in the UHECR 
spectrum seen at Earth. 

The top-down prediction of photon preponderance in 
UHECRs is the one major problem in an otherwise simple 
explanation. Results from the Fly's Eye 0], Haverah 
Park [l^, and AGASA [l4| cosmic ray experiments all 
indicate that at energies above around 10 19 to 10 20 eV, 
protons are more abundant than photons in UHECRs. 
While photons are disfavored, it is not possible to be 
sure that the primary is indeed a proton. 

The idea of violating LI (see Ref. ^| for a broad 
overview) has recently been studied in the context of 
UHECRs for proton decay 0, EJ, atmospheric shower 
development [18| , and modifications to the cutoffs of pro- 
ton and photon spectra due to cosmic background fields 
d ES SJ El El E3- As shown in Ref. [H, LI vi- 
olation can also lead to vacuum photon decay, and in 
general the decay of particles to more massive species. 
We exploit this fact to show how particles produced in 
top-down decays can undergo an "inverse cascade" to 



produce superheavy UHECR primaries. 



II. BREAKING THE LAW 
A. Modified Dispersion Relations 

There are two approaches to brea king LI commonly 
used in the literature. The first E3, El E! El is 
generically to modify the dispersion relations for a parti- 
cle of mass m and 3- momentum p = |p|. 

The second method, which we will not utilize here, is 
instead to write down a particle Lagrangian 27] which 
includes Lorentz and CPT violating terms. A mapping 
of the Lagr ang ian terms to a dispersion relation is non- 
trivial (see |23j for a simple case); however, to first order 
the changes to the photon and electron propagators in- 
duce shifts to c 7 and c e |17|. Similarly, it was shown by 
[28| that loop quantum gravity effects produce modified 
dispersion relations similar to those considered elsewhere 
and here. 

The simplest way to break LI, as shown in Tt\ . is to 
write down a dispersion relation for a particle species i 
as 



E 2 



2 2 i 2 4 

Pi c i + mi c, . 



(1) 



This is changing the "speed of light" or Maximum At- 
tainable Velocity (MAV) for each particle to something 
slightly different than c. 1 The MAVs for different par- 
ticle species are assumed a priori to be different in that 
and similar treatments |26|. 

In such a case, it is possible for previously unallowable 
reactions to occur, such as p — > ne + v e [laUJl' 7 ~~ * e + e~ 



El El and 7 — > vv [25j. A conclusion reached, us- 
ing the first reaction, is that neutrons may make up the 
dominant baryonic component of the UHECR primaries. 



1 Although factors of c are explicitly included in equations, we still 
take it to be dimensionlcss with = c = 1. 
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This, strictly, would be true if one was only tuning the 
proton and neutron MAVs and leaving all others equiv- 
alent (= c). On the other hand, keeping every MAV as 
individually tunable gives an overwhelming number of 
free parameters. 

In order to examine the consequences of varying all 
particle MAVs using the fewest number of free parame- 
ters, we extend the method used in Ref. [2l| and write 
the dispersion relation as 



E 2 = p 2 c 2 + m 2 c 4 + p 2 c 2 f (p/Mc) 

+m 2 c A g (p/Mc) + mc 2 pch (p/Mc) 



(2) 



Here, f(x),g(x), and h(x) are dimensionless universal 
functions having the property f(Q) = g(Q) = h(0) = 0, 
so that as p/M — > the normal dispersion relations are 
recovered. Here M is the mass scale that determines the 
relative degree of Lorentz violation; for our purposes we 
set it at the Planck mass M ss 10 19 GeV. Expanding 
these functions in Taylor series about x = and keeping 
terms of 0(p 2 ) and smaller gives 



B. Example: Photon Decay 

In this section we consider the tree-level photon decay 
process 7 — > e + e~, which is kinematically forbidden for 
c e = c. Following Ref. |17| . we define the parameter 5 ie 



c 



cl 



(9) 

Using astrophysical constraints [2^, , &ye can be lim- 
ited to the range 10~ 16 > 8 ie > -lCT 17 . For 8 ie > 0, 
photon decay occurs above a photon energy threshold 
given by 



E, 



th 



2m e c e 



(10) 



This decay rate of was computed in Coleman & Glashow 
(1997) as 



E 



1 - {E th /Ey 



3/2 



(11) 



E 2 



m c + mc vc —a (0) — 
2 M 



+p 2 c 2 [ 1 



2 K ' M 



1 »/ x m 



(3) 



The g'(0) term can be neglected for two reasons: at 
p <C mc it would be experimentally detected as devia- 
tions in the non-relativistic kinetic energy; for p 3> mc 
it is negligible to the quadratic p term. Note that to 
this (lowest) order, f(x) has no effect. 2 For this section, 
we examine the cases of massless (photon) and massive 
(electron) particles: 



1 



1 



E 2 e = m 2 c*+pjc 2 { 1 + -h'(0)- + ~<7"(0)^ ) (4) 



M (i 



E: 



2 2 



- 7 

We can now define 



(5) 



(7) 



to then write the dispersion relation in the familiar form 
of Eqn. \I\&s 



For decay product energy above about 2x 10 20 eV, pho- 
tons will outnumber protons, so we set E t h to this value 
and suggestively rewrite the decay rate above threshold 



r ( e 

Yc * (1 kmrl (2 x 10 2 ° eV 



3 x 10" 



(12) 



This would rule out all but a terrestrial source of UHE 
photons. The value S ie — 3 x 10 -29 is used is a lower 
limit; if 5 ie would be any lower the energy threshold 
would become too high. This means that only d je G (3x 
10 -29 ,10~ 16 ) will correct for the photon overabundance 
in top-down scenarios. 

Now we see what is required to have 5 ie in the neces- 
sary range to allow for this photon decay. Using Eqs. © 
and J5J|, we can write 



J ye 



1 9, . m> 
-c 2 h'<0) — 

2 V ' M 



1 9 /, / x m 



(13) 



Since the functions g{x) and h(x) are a priori ar- 
bitrary, so are the values and signs of their deriva- 
tives. It is argued in Ref. jH| that too strongly vary- 
ing functions would be unphysical, and that the deriva- 
tives should be of order unity. 4 Adopting this, we can 
ignore the g"(0) term in the above equation. In this 
case, in order to have photon decay we need to require 



1^2 2 4,2 2 

E e = m e c e +p e c e . 



(8) 



2 Issues o f photon dispersion in the vacuum, as discussed in Refs. 
USUI, will not be addressed here. 



3 This is identical to £ — rj for n = 2 in the notation of Ref. I'ii.l . 

4 Note that if this holds also for f(x), we erred in not including the 
/'(0) term in Eq. J^li as that dominates over the h'(0) term when 
P ^ mc - We will not apologize much for this, as we are looking at 
the "lowest order" change in phenomenology (the change in the 
"speeds of light" ) and not "next order" changes (such as vacuum 
photon dispersion). 
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h'(0) negative. Then 5^ 



E, 



2m e 



iy/2Mj~r 



r,2 - 



mc 2 /2M ~ 
2yJ2m e Mc 2 



3 x lfr 23 and 
w 2 x 10 17 eV. 



This value of S^ e is (logarithmically) in the middle of our 
allowable range. 

Returning to our original impetus of addressing the 
feasibility of "Top Down" scenarios, we see that requiring 
h'(0) ~ — 1 is sufficient to eliminate any photons among 
the UHECR primaries. This model, however, eliminates 
any protons as well, which we will now address. 



C. Making the impossible possible 



Unlike previous investigations [13, |29j where the MAV 
for each particle could be chosen independently, this is 
not so in our case. Since Eq. [3 holds for all types of 
particles, it is straightforward to see other effects of this 
level of Lorentz invariance violation. Using Eq.|Bl we can 
in general write for two particles A and X 



6 A 



x 



2 m A - m x 
M 



(mx ~ m A )c 2 
2M 



(14) 



where we have used h'(0) ~ — 1. 

The consequences of this for particle decays are as fol- 
lows, as first noted by 01 > and we reproduce their argu- 
ment here. 

Consider the decay A — > {X}, where {X} is a col- 
lection of massive particles. We want to find the min- 
imum possible energy E min where this can occur. We 
can always lower final state energy by removing trans- 
verse components of momenta, so we're working in one 
dimension only. We want to minimize 



X 



subject to the constraint 



G({px})= (5>x) -PA = 0. 



(15) 



(16) 



x 



The {px} are variables, and pa, Ea > E min are con- 
stants. Using the method of Langrange multipliers, E is 
minimized when 



which becomes 



dp 



Ex 



dE 8G 
— A - 



dp 



x 



Pxc\ 



m x c x+Px c x 



(17) 



(18) 



Thus all products move with the same velocity A. Note 
that A < cx for all X. For the example we consider, in 
the event of massless particles among the decay prod- 
ucts, they would minimize final state energy by carrying 



exactly zero momenta Since the {cx} are the max- 
imum attainable velocies for the decay products, that 
implies that A < minjcjf}. 

It is this last point we focus on: min{cx} in our case 
corresponds to maxfmx}- In the ultrarelativistic limit, 
Ea ~ PaCa and so on, and the particle with the lowest 
speed of light will carry most of the momentum (in the 
minimum energy case). Further, in this threshold case, 
A is also the velocity of the incident particle 0. Since 
A w ca in the regime we are considering, that means that 
the decay can only occur if ca > min{cx}- From Ea. 1141 
that implies that ttia < max{mx}- 

In other words, above a certain threshold, a decay may 
only be allowed if at least one decay product is more 
massive than the decaying particle. This can have serious 
consequences for the particle content of UHECRs, as we 
will now address. 



III. THE INVERSE CASCADE 

In this model, for each individual reaction there is a 
threshold energy Et above which particles can only decay 
if at least one of the decay products is more massive. For 
the decay A — > B + where B is the most massive 
decay product, raj > ttia, and m... <C tub the threshold 
energy goes like 



Using equation 1141 this reduces to 



Eq 



c 2 ^M(m A + m B ) 



(19) 



(20) 



Due to the large value of the Planck mass, this thresh- 
old is around the same order of magnitude for a wide 
range of reactions involving particles of mass much less 
than M. 

Consider an UHECR of mass mo produced with energy 
Eq ^> Et- If normally (at energies E <C Et), particle 
can be a decay product of a more massive particle 1, 
then here (E 3> Et), can decay into 1. We ignore the 
lighter products that will result as well; the substantial 
fraction of the incident energy will be carried by the most 
massive decay product. 

Now, we have a particle of energy E\ < Eq and mass 
mi > mo. If it can, it will also decay into a particle 
2 with m,2 > mi and energy Ei < E\. As the decays 
continue into more and more massive final states, the 
final energy continues to decrease, until at some point 
we reach a particle N with energy En < Et and mass 
mjv > ttiq. Since we are now below threshold, the N 
particle can decay as usual into N-l and so on into less 
massive particles. 

So, to summarize, this inverse cascade decay process 
occurs in two stages. 
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1. Decay up: Decays into more and more massive, 
yet less energetic particles, when E > E T , and 

2. Decay down: Regular decay schemes once E < 
Ex . 



IV. PANDORA'S BOX (NEWLY ALLOWED 
REACTIONS) 

We saw in Section II that at sufficiently high energies 
particles can only decay if one of the decay products is 
more massive. In Section III we saw how this can produce 
an "inverse cascade" to more and more massive particle 
species. In this section we address the specifics of how 
this cascade occurs, by examining which previously for- 
bidden decays can now occur. All of the following are 
assumed to be occuring well above threshold, and we list 
them in several classes. 



A. Tree-level QED 

Here the bare QED vertex 7 — ► Q + Q~ where Q is 
any charged particle is allowed as a decay This also 
includes final states with extra photons and bound states 
of Q+Q-, likc7^7r°7. 



than the top quark. If it is lighter than the top, the Higgs 
is unstable to H — > tt. 



C. Tree-level Strong 

Since the bare QCD vertex preserves quark flavor and 
charge, not much happens in the bare quark sector, apart 
from the reversings of strong decays such as N — > Air and 
7r — ► p-K. Just as for the massless photon, free gluon decay 
is allowed via g — > qq. 

D. One-loop processes 

One could formulate many more decays at the 1-loop 
level, but the rates for these for them are suppressed 
compared to the tree level decays. For this reason, we 
do not consider them further here, except to note two 
interesting decays of the photon: 

7 — ► vivi (QED and weak process). Note that since 
neutrinos are so much lighter than their heavy partners 
and the quarks, the threshold for photon to neutrino de- 
cay could be very low depending on how light the neu- 
trinos are. 

7 — > Z7 (weak process) . This photon Z-Cerenkov de- 
cay occurs due to a W-loop. 



B. Tree-level weak 

It is in the weak sector that most of the interesting 
decays in the inverse cascade occur, due to the flavor- 
changing W vertex. 

In this Lorentz-violating scheme, weak decays (such 
as that of the neutron) can happen in reverse, allowing 
proton decay via p — > ne + v and flavor changing decays 
such astn A-7T . 

Not only that, decays to weak bosons become permis- 
sible, such as vi — > Z ± W T and Z ± — > v\W^ . Included 
in this are decays of leptons to heavier leptons due to 
virtual W's. In the quark sector, this allows decays such 
as p — > nW + and so on. Generally the flavor changing 
weak decay q —* q'W will convert all quarks to t quarks. 
In the absence of a fourth generation, a bare t quark is 
stable. 

Since the Z is massive, decays such as {q,l,fi} — > 
Z{q,l,vi} ("Z-Cerenkov") become permissible as well. 

For the bosons; W + — > tb is allowed, and Z-Cerenkov 
is possible for the W as well (this is a bare weak vertex, 
as is W -> WZZ and W -> WZj). The Z is unstable 
to the decay Z — ► tt. 

Now consider the Higgs: indeed we might be above the 
energy of weak symmetry breaking. Including a single 
massive Higgs H would allow; {/, W, Z} -> {/, W, Z}H 
("Higgs-Cercnkov" for fermions and weak bosons) and 
{W, Z} — > {W, Z}HH. In the absence of new physics 
beyond the Standard Model, the H is stable if it is heavier 



E. Spin up 

For bound states (and now restricting ourselves to 
bound states of three valence quarks) there exist higher 
mass resonances that only differ in the mass and total 
angular momentum (e.g., the A and higher resonances 
for the proton and neutron). Since the state with higher 
angular momentum is more energetic, it is more massive 
and thus processes like N — > N*ir, N* — > N**ir and so 
on can act to "spin up" the particle to more and more 
massive states. 

Thus, even the all top quark baryon ttt isn't completely 
stable, but can decay via spin up. 

V. ESTIMATES OF DECAY RATES 

For the photon decay 7 — ► e + e _ , the decay rate T well 
above threshold was 

r 7 ^ ee = ^aS, /e E. (21) 

We adapt this form of the decay rate to other processes 
to make estimates, and argue for it as follows. At ultra- 
relativistic energies, the only energy scale is the incident 
energy E, so T cx E. It is proportional also to the ef- 
fective coupling constant involved in the decay, and to 
kinematic factors. We schematically write this as 

r <~ {coupling) x {kinematics) x E. (22) 
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For the photon decay, the coupling is a and the kinematic 
factor is 5~ e . We generalize these factors as follows. 

We (first) consider decays of the form A — > B + C 
where ttla < ttib, mB 3> m c- At threshold and beyond, 
most of the momentum of the final state is carried by 
the particle with the lowest speed of light, thus highest 
mass (particle B). The relevant kinematical factor will 
then be Sab- This will continue to be true when final 
states with more than 2 particles are allowed, as long as 
B continues to be the most massive particle. 

The coupling factor is determined by looking at the 
tree level Feynman diagrams and counting vertices in- 
volved in the reaction. 

Consider then the decay of the proton; in addition to 
the inverse of the neutron decay p — ► ne + v e ("leptonic" 
decay), the decay p — > nW + ("bosonic" decay) is also 
allowed. 

For the bosonic decay; Et ~ \pnvwM ~ 2.8 x 10 19 
eV, and as most of the energy goes into the W, S p w ~ 
mw/2M ~ 4x 10~ 18 . The vertex is a single weak vertex; 
with coupling aw = a/ sin 2 9 W ~ 0.0316. Thus, at E = 
£ T ,r~4eV. 

For the leptonic decay; E T ~ • V /M(2m p ) ~ 4.5 x 10 18 
eV, and as most of the energy goes into the n, 5 pn ~ 
(m„ — m p )/2M ~ 5 x 10 ~ 23 . This reaction has two 
weak vertices with coupling a w ~ 10 ~ 3 . Thus, at 
E = Et , r ~ 2.2 x 10" 7 eV. 

Note that the distance travelled before decay is ~ 1/r, 
and given tic = 6.57 x 10~ 24 eV pc, the 1/r for these re- 
actions is incredibly small: 10 -24 pc for the bosonic and 
10 -17 pc for the leptonic. In the context of a UHECR, 
assuming comparable rates (within a few orders of mag- 
nitude) to the above, this inverse cascade happens essen- 
tially immediately, so that by the time any decay prod- 
ucts reach the Earth, they have gone through both stages 
of the cascade and are now sub-threshold particles. 

VI. CONSEQUENCES FOR UHECR 
COMPOSITION 

Given the rapid rate of the inverse cascade process, 
there can be two outcomes for the UHECR population 
at earth. 

First, if both the Decay Up and Decay Down sequences 
occur, then the consituent particles could, depending on 
the exact decay chains and nature of initial UHE primary, 
consist of protons, electrons, photons, and their antipar- 
ticles. Each particle spectrum would exhibit a cutoff at 
their prospective threshold energies, which are all related 
by the h! (0) parameter. 



Second, and more interesting, if the initial UHE pri- 
mary is of high enough energy, the decay up sequence 
may occur and leave particles that can no longer decay 
into more massive particles, but are still above threshold 
for decay into less massive particles. These Most Mas- 
sive Particles (MMPs) would then be the most energetic 
component of the UHECR population. What exactly the 
MMPs are depends on the model of particle physics used. 
Of particles known to exist know, the MMPs would con- 
sist of top quarks and excited bound states of top {ti 
mesons and ttt baryons) quarks. If the Higgs is more 
massive than the top, then the MMPs would consist of 
single Higgs bosons. Including Supersymmetry and its 
expected spectrum of more massive superpartners, if the 
decay rate is still rapid enough the MMP would be the 
most massive superparticle (MMSP). 

Of course, in the first outcome, our proposed solution 
of using LI violation to solve the photon abundance is- 
sues in top-down models is not obviously successful. For 
both outcomes, in fact, one would have to re-examine the 
decay process starting with the initial decay of the super- 
massive relic (which is assumed to be sub-threshold, so 
it decays "normally" ) , and follow the decay products to 
see what results. 

Simulations of MMP-induced shower development 
would have to be performed to see if they are consis- 
tent with actual events. It was shown in Ref. [32 ] that 
stable UHE primaries of mass greater than about 50 GeV 
should be distinguishable from proton primaries by their 
atmospheric shower profiles. As the lightest conceivable 
MMP (a ti meson) would be much heavier than this 
bound, the result of [13] might be used to exclude the 
MMP as the primary for UHE showers. However, a key 
assumption is that the primary is stable. For a MMP just 
above the threshold energy, it may take only a few colli- 
sions before the MMP energy is lowered below threshold, 
where it will decay "normally" (e.g., ti — ► W + W^~) into 
a shower of particles, perhaps mimicking the shower of 
a nucleus UHECR. A MMP with incident energy farther 
above threshold would then penetrate deeper. In any 
event, further investigation into MMP UHECR shower 
profiles is warranted. 

Also of concern is that, while energy and momentum 
are conserved in a certain preferred frame, they are nec- 
essarily not conserved in all other frames. As such, due 
to the motion of the earth and the solar system, small 
departures from energy conservation might be observed 
as well in the UHECR atmospheric showers. 

This work was supported in part by the Department of 
Energy and by NASA (NAG5-10842) and by NSF grant 
PHY 00-79251. 
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